Circuit QED and sudden phase switching in a superconducting qubit array 
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Superconducting qubits connected in an array can form quantum many-body systems such as the 
quantum Ising model. By coupUng the qubits to a superconducting resonator, the combined system 
forms a circuit QED system. Here, we study the nonhnear behavior in the many-body state of the 
qubit array using a semiclassical approach. We show that sudden switchings as well as a bistable 
regime between the ferromagnetic phase and the paramagnetic phase can be observed in the qubit 
array. A superconducting circuit to implement this system is presented with realistic parameters . 

PACS numbers: 



Circuit quantum electrodynamics (QED) has been in- 
tensively explored in the past few years []], 'l] . A number 
of quantum optical effects such as Lamb shift, nonlinear 
spectrum, and lasing have been demonstrated in experi- 
ments [3]. Most recently, the Tavis-Cummings model has 
also been tested with multiple transmon qubits coupling 
to a resonator 0| . The superconducting resonators hence 
not only provide a powerful tool to implement quantum 
information protocols [1-01 , but can also facilitate the 
study of quantum coherence effects and many-body ef- 
fects involving microwave photons jl, Q . 

When a quantum many-body system is coupled to a 
resonator cavity, novel many-body phenomena that have 
not been studied previously in condensed matter systems 
can be observed. In the atomic systems flO], it has been 
shown that atomic condensate can be coupled strongly to 
cavity photons to demonstrate interesting effects (ijj, |l2| . 
In the superconducting circuits, given the diversity of 
the electromagnetic couplings and parameter regimes, a 
large variety of many-body Hamiltonians can be engi- 
neered that can not be emulated in atomic systems. In 
the past, critical behavior has been studied both theoret- 
ically and experimentally in Josephson junction arrays 
that can emulate spin models and quantum phase mod- 
els [l3|- Recent progresses in superconducting quantum 
devices further expand such research to include fermionic 
Hamiltonians [14]. 

In this work, we study the nonlinear behavior in the 
many-body state of an array of superconducting qubits 
coupled to a superconducting resonator. Different from 
previous work on optical bistability [15, 16], our work fo- 
cuses on the behavior of the qubit array in this system. 
The qubits are arranged in a one-dimensional chain to 
form the quantum Ising model. In the bad cavity limit, 
the photon number of the resonator will depend nonlin- 
ear ly on the average of the qubit operator. This gives 
rise to a sudden switching between different quantum 
many-body phases in the qubit array when the driving 
power on the resonator is varied, instead of the continu- 
ous quantum phase transition in a simple quantum Ising 
model. Our study also shows that a bistable regime exists 
where the qubit array can be in either the paramagnetic 



phase or the ferromagnetic phase. A circuit to implement 
this model can be constructed with realistic parameters, 
where the inductive coupling between the superconduct- 
ing qubits and the resonator generates a shift of the res- 
onator frequency that depends on the qubit operators. 
The nonlinear effects studied here can be experimentally 
demonstrated in a toy system with only a few qubits and 
a resonator. Furthermore, the resonator can act as a 
threshold detector to measure the phase switching in the 
bistable regime [l3| and act as a probe of the continuous 
quantum phase transition in the weak driving limit [18]. 

We consider a qubit array coupled to a resonator with 
the total Hamiltonian 



Ht = Hqi + H^ 



int 



H^ 



cO 



H^ + Hy 



(1) 



where Hqj = -JxY^i^xi - Y^i^ziCTzi+i describes the 
quantum Ising model with a uniform transverse mag- 
netic field Jx and a ferromagnetic coupling between near- 
est neighbor qubits, i^co = huca^a + he{t){a + a"*") de- 
scribes the Hamiltonian of the resonator and includes a 
driving term with the amplitude e{t) = 2ecoscjt, and 
Hint = ^gcL^CL^cFxi is the coupling between the res- 
onator and the qubits. Here, Gxi and cFzi are the Pauli 
matrices of the qubits and a [a^) is the annihilation (cre- 
ation) operator for the resonator. The magnitude of the 
ferromagnetic coupling is set to unity and is used as the 
energy unit in the following. The term H^ {H^) is the 
coupling between the resonator (qubit) and its environ- 
mental modes that causes damping (relaxation) of the 
resonator (qubit). Without the resonator, a second or- 
der quantum phase transition between the paramagnetic 
and the ferromagnetic phases occurs in the qubit array 
when the transverse field is varied [l9]. 

The coupling Hint shifts the frequency of the resonator 
by ^ ^ axi. In the Heisenberg picture, we can derive the 
operator equation for the resonator in the rotating frame. 
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where Ac is the detuning, n is the damping rate of the 
resonator, and a^n is the noise operator (l5|. In the bad 
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FIG. 1: (a) photon number Us vs. driving power e . Solid 
curves from left to right are for Jx = 1.4, 1.6, 1.8, 2.0 at Ac — 
0; dotted curve is for Jx = 1.8 at Ac = gM. (b) X vs. for 
Jx = 1.8 at Ac = (solid curve) and at Ac = gM (dotted 
curve); inset is X vs. Jx in the quantum Ising model, (c) 
stability coefficient Cs vs. for Jx = 1.8 at Ac = 0. We use 
K = 0.03, and g = 0.0005. The dashed lines indicate ei^2- 



cavity limit with strong resonator damping, we can derive 
the operator average in the steady state as 



(a) 



{iAc - n/2) - ig{Y.Mxi))s 



(3) 



using a semiclassical approximation with (a ^i{cFxi))ss ^ 
{'^) ss(^i{^xi)) SSI where the quantum correlation be- 
tween the qubits and the resonator is neglected [16]. The 
photon number in the steady state is then 



Us = —{a — o/) 



t\ - 



(4) 



depending nonlinear ly on (X]i(^a^i))ss- considering 
the relaxation of the qubit array, we can take the steady 
state average to be {'^i{o'xi))ss = ^{Jx), where X is the 
ground state average of ^^{crxi) in the quantum Ising 
model with the effective transverse field Jx = Jx — ^g^s 
and can be calculated with the Jordan- Wigner transfor- 
mation. Due to the back-action of the resonator, the 
effective transverse field is now shifted by a term propor- 
tional to the photon number. 

Below we will study how the nonlinear dependence af- 
fects the resonator and the quantum many-body system 
[20]. Choosing Ac = and J^^ > 1 to be above the crit- 
ical point of the quantum Ising model Jc = 1, we first 
plot the photon number Ug versus the driving power 
in Fig. [TJa). At low driving power, Ug is small so that 
Jx ^ Jx > ^ and the qubit array is in the paramagnetic 
phase with X ^ M, where M is the number of sites in the 
array. The photon number increases nearly linearly with 
e^. As the driving is increased, becomes larger and Jx 
decreases to approach the critical point. In this regime, a 
bistable behavior can be observed: when is increased. 
Us initially follows the lower branch but suddenly jumps 
to the upper branch at e^; when is decreased from an 
initial high driving power, ris follows the upper branch 



but suddenly jumps to the lower branch at . Mean- 
while, the jump in the photon number is accompanied 
by a sudden switching of the quantum many-body state 
in the qubit array between the ferromagnetic phase and 
the paramagnetic phase which can be seen in Fig. [TJb). 
Hence, in the bistable regime between ef and 63, the qubit 
array can be in either the paramagnetic phase or the 
ferromagnetic phase. Note that the bistability depends 
critically on the system parameters. For a detuning of 
Ac = ^M, no bistable regime can be found. 

The stability of the semiclassical solutions can be de- 
termined from the dynamic equation for the photon num- 
ber: dn/dt = ie{a — a^) — K,n [2u\. Assume a small fluctu- 
ation Sn in the photon number n = ns-\-Sn. The average 
X obtains a fluctuation 5X = —hgX'Sn with X' being 
the derivative of X. And the fluctuation S{a — a^) can 
be derived from Eq. (|3|). We can then derive 



dSn 
~dt 



— — = —ndnil 



2nl^hg" 



X'X) 



(5) 



at Ac = 0. When the stability coefficient satisfies Cs = 
1 — 2nl^hg^X' X/ > 0, the semiclassical solutions are 
stable [20]. In Fig. [IJc), Cs is plotted which shows that 
the solutions in both the upper branch and the lower 
branch are stable while these in the middle branch are 
unstable with Cg < 0. _ 

The effective field J^^of the qubit array also makes a 
jump from J^^ > 1 to J^; < 1 at 62 and vice versa at ef 
as is shown in Fig.[2ja), when the photon number makes 
the sudden jump. This shows that the qubit array under- 
goes a sudden switching between the paramagnetic phase 
and ferromagnetic phase at 2? and does not cross the 
point Jx = ^ continuously. Meanwhile, the ground state 
energy of the qubit array obtains a finite difference as 
well. In Fig.[2](b), we plot the energy difference per qubit 
AEg/M at the switching points 2 respectively, where 

the ground state energy of the qubit array at the field Jx 
is derived using the Jordan- Wigner transformation: 



(6) 



with E{x) being the complete elliptic integral of the sec- 
ond kind. Such behavior in this driven system is in anal- 
ogy to the first order phase transition in an equilibrium 
system. While in a simple quantum Ising model with- 
out the resonator, continuous quantum phase transition 
occurs at the critical point Jc = 1. 

The boundaries of the phase switchings 2 ^ func- 
tion of Jx can be derive by seeking local maxima and local 
minima in the curve of versus for each Jx. Using 
the phase boundaries, we can characterize the quantum 
many-body phases of the qubit array in terms of two pa- 
rameters: the transverse field Jx and the driving power 
e^, as is shown in Figl2](c). For > 63 (J^^), the system 
is always in the ferromagnetic phase; for < (J^^), the 
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FIG. 3: Circuit of an array of superconducting qubits coupling 
to a resonator. The current in the qubit loops couples to the 
SQUID in the resonator via mutual inductance. 



FIG. 2: (a) Jx vs. Jx before the sudden switching (solid curve) 
and after the switching (dashed curve) at e\^2- (b) the energy 
change at the switching points vs. Jx- (c) many-body phases 
in the space of Jx and e^, where 'F' (horizon lines) indicates 
the ferromagnetic phase and 'P' (vertical lines) indicates the 
paramagnetic phase. 



system is always in the paramagnetic phase. In the over- 
lapping area between these two boundaries, both phases 
can exist depending on the history of the parameters. 
The coupled system of the qubit array and the resonator 
can thus be explored to demonstrate nonlinear effects in 
the quantum many-body state. 

The quantum Ising model can be constructed with su- 
perconducting qubits [Hll connected in a chain and cou- 
pled capacitively to their neighboring qubits, as is shown 
in Fig. [3l The Hamiltonian can be written as 



Hqi 



'^Eji cos{ipi) 



(7) 



where = [rii — {CgiVgi/2e), ■ 



is the number of Cooper pairs on the superconducting 
islands biased by the gate charge number CgiVgi/2e^ and 
is the inverse of the capacitance matrix with Cu = 
Co being the total capacitance connected to the i-th is- 
land and Cii±i = —Ci being the coupling capacitance 
between two adjacent islands with Ci < Cq. Expanded 
to the second order of Ci/Co, the elements of can 
be derived as: C^^ = I/Cq + 2Cf/Cl 



and C- 



i+n 



o{Ci/C^){l/Co) for n > 1. To obtain the 
ferromagnetic coupling, we define the spin states of the 
i-th qubit by the charge states |0)i and on the i-th 
island as: |0)i \ t)*,]!)* ^ | i)i when i is an odd 
number; and ^ | t)^, \0)i \ Di when i is an even 
number. By setting the gate voltages to be at the degen- 
eracy point with CgiVgi/2e = 1/2, the single qubit term 
in Eq. ([7j) becomes '^{—Eji/2)(Jxi'-, and the coupling 
term can be written as ^^(-^licr^icr^i+i + B2icrzicrzi^2) 
to the second order of Ci/Cq with Bu = 4e^(Ci/Co) and 
B2i = 4e^(Ci/Co)^. Due to the long range nature of the 
Coulomb interaction, the coupling term includes not only 
the nearest neighbor couplings, but also unwanted cou- 
plings between other qubits. However, it has been shown 
that the phase diagram of this model is qualitatively the 
same as that of the simple quantum Ising model when 



B2i/Bu = Ci/Co < 1/2 [27|. In our system, Ci/Cq can 
be readily designed to satisfy this criterion. 

To generate a coupling Hint = hga^a^cjxi between 
the qubits and the resonator, we consider a resonator 
that includes a SQUID in the circuit. In Fig. [3l each 
qubit has a small loop current Iqi = Iq2 ^\ii{52/2)crxi 0| 
that is controlled by the external flux 82- This cur- 
rent is coupled to the SQUID via mutual inductance 
and modifies the flux in the SQUID loop by = 
— /io^g2^o sin((52/2) (Ja^^ where Rq is roughly the size 
of the qubit loop. The effective inductance of the SQUID 
can then be written as 



(-) 



1 



2e 2Ir cos (^ex — 2Ir siu (pex {TT^r/^o) 



(8) 



where Ir is the critical current of the SQUID junctions 
and cpex is the external flux in the SQUID loop. The total 
inductance of the resonator is Lfot = + ^r, including 
the SQUID inductance and the loop inductance L^. Let 
L^q be the SQUID inductance at <l>^ = and cp^x = 7r/4, 
the frequency of the resonator becomes 



^sq 



(9) 



with Uco = {{L^q -\- Lr)Cr)~^^'^ ' This expression includes 
the coupling between the qubits and the resonator in the 
form of Hint with 



= ^co( 



^sq 



Lr 2 ^oV2 



■)• 



(10) 



The coupling constant g depends on the qubit flux to 
the first order, which is a merit of this circuit that ensures 
a sizable coupling to demonstrate the nonlinear effect. 
We can choose the following parameters: Lr ~ 100 pH, 
Iq2 = 80 nA, Ir = 1200 nA (i.e. L^^ = 100 pH), and the 
total capacitance of the resonator circuit Cr = 0.1 pF, 
which give ^ = 27r x 1 MHz and Uc = 27t x 29 GHz. With 
a ferromagnetic coupling of 27r x 2 GHz chosen as the 
energy unit, we have g = 0.0005 in the relative unit. Note 
that various other qubits such as the superconducting 
flux qubit and phase qubit can be used to construct the 
qubit array given the flexibility in the superconducting 
circuits [23.-.25J. 
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Decoherence plays an essential role in observing the 
nonlinear effect. In the above, we assume that the res- 
onator is in the bad cavity limit where the damping of the 
resonator is stronger than the coupling constant and the 
decoherence of the qubits. This can be achieved either 
by using a noisy resonator or by connecting the resonator 
to a dissipative element which is shown as the resistor Rr 
in Fig. [3l The decoherence of the qubits will relax the 
quantum many-body system to its ground state at the 
effective transverse field. Given the sub-megahertz de- 
coherence rates that are recently measured in the super- 
conducting qubits, the nonlinear effect can be studied in 
the above circuit with realistic parameters. Finite tem- 
perature may play a role in this system as well when the 
field Jx becomes much smaller than the ferromagnetic 
coupling, where a low-lying excited state is nearly degen- 
erate with the ground state. The final state can then be 
a mixture of the ground state and the excited state at 
a temperature of T = 20 mK {kBT/h = 27r x 0.4 GHz). 
In addition, the detection of the phase switchings in the 
qubit array can be achieved by measuring the photon 
number in the resonator. The jump of the photon num- 
ber is directly associated with the phase switching when 
the driving is in the bistable regime. 

To conclude, we studied the nonlinear effect in the 
quantum many-body state of a superconducting qubit ar- 
ray coupling to a superconducting resonator. We showed 
that sudden phase switching between the paramagnetic 
phase and the ferromagnetic phase can occur in the qubit 
array when the driving on the resonator is varied due to 
the strong back-action of the resonator on the qubits. 
The phase switching results in a finite energy jump in 
the qubit array, in analogy to a first order phase transi- 
tion. These results showed that novel phenomena can be 
observed when a global resonator cavity is connected to 
a quantum many-body system. Many interesting effects 
await to be discovered in such systems. 
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